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ABSTRACT 

Using the results of our first paper on the Chandra HRC observation of the Orion 
Nebula Cluster (ONC), here we explore the relation between the coronal activity of its 
1-Myr-old pre- main sequence population and stellar parameters. We find that median 
X-ray luminosities of low mass stars (M/Mq < 3) increase with increasing mass and 
decreasing stellar age. Brown dwarfs (0.03 < M/Mq < 0.08) follow the same trend with 
mass. From M ~ 0.1 to M ~ 0.5 M@, median Lx/Pboi values increase by about half an 
order of magnitude and then remain constant at ~ 10~ 3 ' 5 for the mass range from 0.5 
to 3.0 M/Mq. In these same two mass ranges, Lx/Lbol remains roughly constant with 
age, until it drops by more than two orders of magnitudes at the epoch when ~ 2 — AMq 
stars are expected to become fully radiative. 

We find a dependence of Lx and Lx/Lboi on circumstellar accretion indicators 
and suggest three possible hypotheses for its origin. In spite of improved X-ray and 
rotational data, correlations between activity indicators and rotation remain elusive 
for these stars, possibly indicating that stars for which rotational periods have been 
measured have reached some saturation level. 

Our study of X-ray activity vs. stellar mass leads us to propose that the few HRC 
X-ray sources not associated with any optical/infrared counterpart trace a yet to be 
discovered stellar population of deeply embedded, relatively massive ONC members. 
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1. Introduction 

The present work focuses on the study of the X-ray activity of Orion Nebula Cluster (ONC) 
members: from a purely observational standpoint, we explore the relationship between X-ray activ- 
ity and stellar characteristics, searching for the physical mechanisms responsible for coronal X-ray 
activity in pre-main sequence (PMS) stars. In spite of the previous observational work in this area 
(e.g., Feigelson et al. 1993; Casanova et al. 1995; Gagne et al. 1995; Flaccomio et al. 2000), no 
definitive picture has yet emerged. 

Our investigation is based on the original Chandra High Resolution Camera (HRC, Murray 
et al. 2000) X-ray data and literature optical data presented in Flaccomio et al. (2002, hereafter 
Paper I) and Flaccomio (2002) for objects in the ONC area. For the ONC region under study, 
Paper I defined an optical sample of 696 optically selected, extinction limited, and well characterized 
ONC members. Here we address the controversial questions of the relationship of X-ray activity 
with convection and rotation (i.e., the classical a — u> dynamo parameters), with stellar mass and 
age, and with circumstellar accretion and/or the presence of circumstellar disks. Taking advantage 
of our detailed description of activity as a function of stellar mass, we then explore the nature of 
the few detected X-ray sources not associated with any optical/infrared (IR) counterparts. 

The structure of this paper is as follows: In §2 we study the relation of activity with rotation, 
mass, age and circumstellar accretion strength. In §3 we speculate on the nature of X-ray sources 
with no optical/IR counterpart. Finally we discuss our findings and summarize our conclusions in 
§4. An Appendix details our technique for investigating possible sources below our X-ray sensitivity 
limit. 

2. Relationships between X-ray activity indicators and stellar parameters 

In exploring relationships between X-ray activity and stellar characteristics through empirical 
correlations with physically meaningful stellar parameters, our main goal is to understand the ~5 
orders of magnitude spread that we observe in Lx and Lx / ^boi ■ We focus on how activity is related 
to stellar rotational period, to mass and evolutionary status (as inferred from evolutionary tracks - 
see Paper I), and to the presence of an accretion disk. Paper I presents relevant X-ray and optical 
data for our optical sample comprised of 696 ONC likely members, 342 of which are unambiguously 
detected in our HRC X-ray data and have both low optical extinction (Ay < 3.0) and available 
mass estimates. This stellar sample was carefully selected in order to avoid X-ray selection effects. 

Our non-parametric approach is to split our reference sample into several bins of the stellar 
parameter under examination, to evaluate for each of these bins a Maximum Likelihood distribution 
function (using the Kaplan-Mayer estimator) for our two X-ray activity indicators, and then to 
examine the trends of these distributions as a function of the parameter. The Kaplan-Mayer 
estimators were computed using the ASURV package (Feigelson & Nelson 1985), which was also 
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used to perform various two population tests 2 on pairs of selected stellar subsamples in order to 
establish the confidence with which we can reject the null hypothesis (that two distributions of a 
given parameter are drawn from the same parent distribution). 

2.1. Rotational Period 

We first searched for a correlation between X-ray activity and stellar rotational period (P ro t)- 
Such a correlation is known to hold for older, main sequence stars of low mass and moderate 
rotation rate and is the main evidence for activity being driven by an a — uj stellar dynamo. Such 
correlations have not previously been found for the ONC stars (Gagne et al. 1995) and several other 
PMS associations (e.g., Walter & Barry 1991; Feigelson et al. 1993; Alcala et al. 1997). In the 
young Taurus region, to the contrary, a correlation between activity and rotation has been observed 
(Bouvier 1990; Strom & Strom 1994; Damiani & Micela 1995; Stelzer & Neuhauser 2001). Given 
the incompleteness of rotational databases used in most previous studies, the question remains 
open. 

In Figure 1 we show the distribution of rotational periods (upper panel) collected in Paper I 
from the recent literature for 282 stars in our optical sample. Also shown vs. P ro t is the HRC 
detection fraction (lower panel) . Our database of rotational periods is significantly more complete 
than that available to Gagne et al. (1995), but in spite of this improvement and Chandra's higher X- 
ray sensitivity (which yields a high detection fraction), we still can find no evidence of correlation 
between X-ray luminosity and P ro t- The scatter plot of Log(Lx/Lboi) vs - Log(P ro t) (Figure 2) 
similarly shows no evidence of a correlation with P ro t- These results hold true no matter what 
physically meaningful subsample of the whole we consider. Indeed, in order to minimize stellar 
structure differences throughout the sample that might wash out a correlation with P rot , we repeated 
the analysis in several mass bins (see §2.2) and also for high and low accretion stars (see §2.5), 
obtaining the same negative results in all cases. 

It should be noted that the sample of stars for which rotational periods are available is strongly 
biased toward the most X-ray active stars. For example, the median Log{Lxl^boi) for stars with 
known rotational period and mass smaller than 3.0M Q is -3.37 while the corresponding median for 
stars of the same mass with unknown rotational periods is -3.73. Moreover, the subsample of stars 
with known P rot appears to be dominated by the most weakly accreting stars in the whole sample 
(see §2.5). 



2 Gehan's generalized Wilcoxon test using both permutation and hypergeometric variance; Logrank test; Peto and 
Peto generalized Wilcoxon test; Peto and Prentice generalized Wilcoxon test. 
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Fig. 1. — Upper panel: Rotational period distribution for the optical sample, the subset of stars 
detected in our HRC data and two subsets of the optical sample segregated by their Ca II line 
equivalent width (see §2.5). Lower panel: detection fraction as a function of P ro t- 



2.2. Stellar mass stars 



Figure 3 shows Maximum Likelihood X-ray luminosity functions (XLFs) for ONC members in 
eight mass bins, while Figure 4 shows a scatter plot of Lx vs. mass with characterizations of the 
XLFs superposed for the same mass bins (see caption). In Figure 4 a direct correlation between 
median Lx and mass can clearly be seen for M < 3M , while in the mass range from 3 to 10 M Q 
(where 6 of 12 stars have luminosities below detection threshold) the median Lx drops significantly 
and only an upper limit can be estimated. For massive O-type stars with M > 10 M , Lx rises 
again. 

Scenarios in which a dynamo drives X-ray activity would interpret this break at ~3 M© as 
the point at which stars of the ONC's age (~ lMyr) become fully radiative and the dynamo stops 
functioning. This interpretation is indeed supported by the stellar structure evolutionary tracks of 
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Fig. 2. — Log(Lx / Lboi) vs. Log{P ro t) scatter plot. Filled circles represent detections, small trian- 
gles upper limits in Log{Lx / L^i). 

Siess, Dufour, & Forestini (2000, hereafter SDF). 

The drop in Lx at ~3 M is even more pronounced than shown in Figures 3 and 4 if we take 
stellar multiplicity into account. The fact that most secondaries will have significantly lower masses 
than their primaries, together with the mass - Lx relation, means that multiplicity effects on the 
XLFs can generally be expected to be rather small. But if primaries in the 3 — 10 Mq bin have very 
little intrinsic X-ray emission, we will detect only secondary emission. Indeed, we observe that the 
XLF for 3 — 10 M mass stars crosses those for stars of lower mass and is determined for about 
half the bin population. This lends indirect support to the hypothesis that emission for 3 — 10 M 
mass stars is mainly due to unresolved secondaries. 

Figure 5 shows distribution functions for Lx / L^oi , in the same mass bins as defined for Figure 
3, while Figure 6 follows the format of Figure 4 to show Lx/Lf, i vs. mass (logarithmically). 
The drop in activity at 3M is even more striking with Lx/L^i as an indicator: more than 2 
orders of magnitude (note that the median value in the 3 — 1OM range is an upper limit). It 
is also interesting to note that the upper envelope of log(Lx / L^oi) for stars with M < 3M is 
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Fig. 3. — Maximum likelihood X-ray luminosity functions of the ONC optical sample, for eight 
ranges of mass from 0.1 to 50 M Q . Note that for some distributions, low luminosity tails are not 
probed due to insufficient sensitivity: the faintest stars in such bins were not detected. 



compatible (within uncertainties) with the level of ~ —3, i.e. the saturation level found for both 
rapidly-rotating main sequence stars (Vilhu & Walter 1987) and low mass stars in several other star 
formation regions (e.g., Alcala et al. 1997; Flaccomio et al. 2000). This is also apparent from Figure 
5. No satisfactory theoretical explanation for the physical nature of this saturation has been found 
to date. Our analysis quantifies typical Lx/L^i values in terms of the median and the scatter of 
points and finds somewhat lower values than the canonical -3 value, indicating either a different 
saturation level or perhaps the presence of a mix of saturated and non-saturated stars. There also 
seems to be a relationship between median level and stellar mass, with the lowest mass stars (e.g., 
those in the O.IO-O.I6M0 bin) having a significantly lower median value of Lx/L^i than those of 
higher mass stars. The range of the variation between the O.1O-O.16M and the 0.5-1.0M Q bins is 
0.45 dex, larger than our estimated uncertainty in individual Lxj^boi values (cf. Paper I). A viable 
theoretical explanation of activity in the PMS phase will have to account for this mass dependence 
of either the fraction of saturated vs. unsaturated stars or of the saturation phenomenon itself. 
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Fig. 4. — Scatter plot of Log(Lx) vs. Log(M/M & ) for the optical sample, with filled circles denoting 
detections and small downward-pointing triangles, upper limits. Typical la and 50% uncertainties 
are given by thick and thin error bars on two illustrative data points at lower right. Also indicated 
are medians (diamonds centered on thick segments) with uncertainties (shaded boxes), and the 
25% and 75% quantiles (thin- lined boxes). Upper limits corresponding to the lowest detection in 
a given mass range are substituted for undefined quantities: a diamond with downward-pointing 
arrow indicates an upper limit to the median; two small downward-pointing arrows at box bottom 
indicate an upper limit on the lower error in median and/or the lower quartile. Numeric values 
above abscissa for each box indicate numbers of detections ("Det.") and upper limits ("Upp.") 
contributing to the XLF calculation. 



To conclude our study of X-ray activity vs. stellar mass, we now wish to address the X-ray 
properties of the significant fraction of undetected very low mass stars. More specifically we wonder 
whether non detected stars are simply located at the bottom of the Lx distributions shown in Figure 
3 or whether they constitute distinct stellar populations with significantly lower X-ray luminosities. 
Appendix A describes how we used our Chandra HRC data to derive (or constrain) mean X-ray 
luminosities for mass segregated groups of objects, some or all of which are individually undetected. 
The method measures the flux of a "composite" source by summing HRC photons detected at the 
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Fig. 5. — Maximum Likelihood distribution functions for Lx/Lf, i, similar to Figure 3. 

optical positions of individually undetected objects comprising the "composite." 

Figure 7 shows the logarithm of average Lx, both for the whole reference sample 3 and for 
undetected objects in several mass bins. The height of the boxes reflect the uncertainty in the 
count-rate to flux conversion factors for the "composite" source (see Appendix A). The number 
of summed objects in each mass bin is reported in the upper part of the figure, and for reference, 
the median Lx from Figure 4 is also shown. As a "sanity check" we compared the average Lx 
computed with the method described in Appendix A for the whole reference sample of detected 
and undetected objects with the average Lx computed from the maximum likelihood luminosity 
functions (not shown in the figure). We find the two results consistent, lending support to this 
analysis. Note that the discrepancy of mean and median Lx in Figure 7 is indeed due to a real 
difference between the two indices and indicates an asymmetry in the Lx distributions. 

From Figure 7 we therefore draw the conclusion that undetected stars in the stellar samples 



3 Note that this reference sample contains fewer objects than the complete optical sample. This is due to a number 
of constraints imposed when selecting objects to which our method was applied (cf. Appendix A) . 




Fig. 6. — Log plot of Lx/Lb i vs. mass, using conventions detailed in caption of Figure 4. 

(0.1 < M/Mq < 0.5) have an average Lx well below our sensitivity limit (note that the 0.1-0.16M Q 
bin has only an upper limit). This could either mean that the XLFs have very long low-luminosity 
tails or indicate the existence of a stellar population with quiescent/suppressed magnetic activity. 

The same qualitative results found in this section are also obtained if the evolutionary tracks 
of D'Antona & Mazzitelli (1997) are substituted for those of SDF in estimating masses. Similarly, 
if we exclude the few optical sample stars with unknown astrometric membership (cf. Paper I) from 
our analysis, the results also remain unchanged. 



2.3. Brown Dwarfs 



Our optical sample contains 15 substellar mass objects (0.008 < M < 0.08 M@) all of which 
were studied through IR spectroscopy by Lucas et al. (2001). Only one of the brown dwarfs (BD), 
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with an estimated mass of M = 0.074 Mq, is detected in our HRC data . Our XLF analysis 
therefore cannot be extended to such low masses. This poses the question of why so few brown 
dwarfs are detected: Are BDs subluminous in X rays with respect to our mass - Lx relation, or 
do they fall on this relation but below our limiting sensitivity level, as a linear extrapolation of the 
mass - Lx relationship in Figure 4 indeed suggests. 

In order to investigate this question, we once again applied the method of Appendix A: the two 
substellar mass bins in Figure 7 show our results. The lowest mass bin contains just 3 undetected 
objects and yields only an upper limit. Nine objects lie in the brown dwarf bin (0.03 < M Q < 0.1) 5 
and our method yields a positive detection of the "composite" source with high confidence 6 . Figure 
7 indicates that: 1) brown dwarfs in the ONC do not depart from the Lx - mass relationship we 
observe for higher mass stars, and 2) the mean BD X-ray luminosity is of the order of 10 28 ' 5 
erg-s -1 . On average this is higher than the recent results obtained for BDs in IC 348 by Preibisch 
& Zinnecker (2001), but lower than the Lx reported by Imanishi et al. (2001) for one detected BD 
and one BD candidate in the young p Ophiuchi cloud. 



2.4. Stellar age 

Although the relationship with mass explains a good deal of the scatter in activity levels of our 
sample, the deviations from this trend are still significant. The next stellar parameter we investigate 
in order to understand this residual scatter is the age, as inferred from the SDF evolutionary tracks. 
Figure 8 shows the Log(Lx) vs. Log(Age) plot for stars in our sample with mass between 0.5 and 
1.0 Mq. The median Lx appears to decrease with increasing age, but we must exercise special 
care in the interpretation of this trend: there is a widespread concern that the age spread indicated 
by the position of low mass stars still on the Hayashi tracks may be largely, or even entirely, due 
to an artificial spread in the bolometric luminosities (cf. Hartmann 2001). The sources of error 
that may contribute to uncertainties in Lb Q i and consequent age are indeed numerous and include 
uncertainties in: adopted distance, spectral type and extinction, photometric variability, unresolved 



4 The brown dwarf status of this object is uncertain. The quoted mass is taken from Lucas et al. (2001). However, 
an independent mass estimate from the optical data of Hillenbrand (1997) and the SDF tracks (cf. Paper I) gives a 
larger value of 0.11 M©. 

5 Although this mass range nominally extends above the BD limit, the most massive object in this bin has a mass 
of only 0.08 M©. 

6 An integrated flux of ~40.1 photons in a total effective exposure time of 553 ksec is measured. As discussed above, 
one of these 9 objects is detected individually, with 16.2 extracted counts. If we exclude this object of uncertain 
BD status, the resulting "composite" source falls slightly below our detection threshold (3a), with its probability of 
being due to background fluctuations at 0.28% (2.8ct). Figure 7 indicates it as an upper limit corresponding to less 
than 25.1 photons accumulated over an effective exposure time of 491 ksec. Had we instead accepted it as a positive 
detection, an accumulated photon count of 23.9 would imply an X-ray luminosity virtually identical to the upper 
limit value. 
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Fig. 7. — Mean and median Lx vs. mass for a subsample of optical sample stars and brown dwarfs 
in different mass bins (see text for details). Hatched boxes (small boxes appear as mere segments) 
represent mean Lx derived for the whole sample (black boxes) and for the subset of undetected 
objects (gray boxes); see Appendix A. Box widths reflect mass bin sizes, while heights reflect 
conversion factor (CF) uncertainties by indicating values obtained with the mean and the median 
CFs of the individual objects within the bin. Stars and upper-limit triangles represent Lx medians 
from Figure 6. Rows of numeric values at top give numbers of objects used to compute averages 
for each mass bin. 
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companions, and accretion luminosity. On the other hand, if one were to concede a real age spread 
in the ONC (plausible as star formation is still taking place at the present day), stellar ages might 
carry a statistical and relative significance, though still uncertain on an individual basis: stars that 
we place lower on the Hayashi tracks may be older, on average, than those that we place higher. 
This latter hypothesis is supported by the time evolution of ONC stellar radii recently observed 
by Rhode et al. (2001). Although our correlation of Lx with age also seems to argue in favor of a 
real age spread, it could also be an artifact of an interrelation between inferred Lx and Lb Q i, and 
the matter remains open. Indeed the trend in Figure 8 corresponds to a constant ratio of X-ray to 
bolometric luminosity at different ages and can be interpreted equally well in two different ways: 
1) Stars of equal mass have equal bolometric and X-ray luminosities, but there are effects that act 
in exactly the same way on both Lx and L^; 2) These stars are saturated and stay saturated 
through their contraction on the Hayashi tracks, thus keeping Lx/L^i constant. 

Figure 9 presents a plot of Log(Lx) vs. Log(Age) for stars in the 2 — 3 M mass bin and 
reveals an age dependence and a sudden drop of Lx at Log(Age) ~ 6.5 (a similar drop is also seen 
in Lx/Lboi). This is just the age at which a 2.5 M & star dissipates its convective envelope according 
to the SDF models, lending further evidence for the convection - activity connection found in §2.2. 
Figure 9 also supports the existence of a real age spread in the ONC. Contrary to the situation in 
the 0.5 — 1.0 Mq mass range, evolution of these higher mass stars occurs in temperature as well as 
in luminosity for Log(Age) > 6.0. The spread in age therefore cannot be due to a spurious spread 
in Lb i, as may be the case in the lower mass bin. 

Figures 10 and 11 show trends of median Lx and Lx/L^oi as a function of mass in each of 
four age bins. Care must be taken in interpreting these plots because many of the mass/age bins 
for which medians are computed are sparsely populated (numbers of detections and upper limits 
are given in Table 1) and some points are subject to rather large uncertainties. In those bins with 
less than four detections, we show means computed from the XLFs (instead of medians), and both 
medians and means are replaced by upper limits in some cases. From left to right in the lower 
portion of Figure 11 we show mass ranges over which stars of each age bin are fully convective, 
partially convective, and fully radiative according to SDF tracks. We can draw several conclusions 
from this presentation: 1) For all ages, Lx increases with mass for low masses and then suffers a 
break at higher mass; 2) The mass at which this break occurs appears to depend on the age, being 
lower for the oldest stars in our sample; 3) There is good agreement between the mass at which the 
break occurs and the mass at which stars in that age bin become fully radiative; 4) The increase 
of Lx/Lboi with mass appears again at the low mass end, at least for Log(Age) < 6.0. 



2.5. Accretion/disk indicators: Ca II 

It has long been disputed whether the presence of circumstellar disks and/or circumstellar 
accretion plays a role in determining the activity levels of PMS stars. We are now in the position 
to address this problem in the ONC. We will use the A = 8542 A equivalent width (EW) of the 
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Fig. 8. — Log(Lx) vs. Log(Age) for stars of mass 0.5-1.0 M ; symbols same as Figure 4. 




Log(Age) 



Fig. 9. — Log(Lx) vs. Log(Age) for stars of mass 2.0-3.0 M ; symbols same as Figure 4. Medians 
for the two left-most bins have been replaced by averages. 
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Fig. 10. — Log(Lx) vs. Log(M/M Q ) in four age ranges. Black- and gray-colored symbols, con- 
nected with different-styled lines, give results for age ranges designated by the legend at upper left. 
Large symbols indicate median Log(Lx) values from Maximum Likelihood distributions in each 
mass range; small symbols indicate mean Log(Lx) for those mass/age bins with < 4 detections. 
Downward-pointing arrows indicate upper limits on medians or means. 

Ca II triplet as an indicator of disk accretion. In main sequence and weak-line T Tauri stars, 
where accretion is absent, this line is seen in absorption with EW ~ 3, a value that depends only 
weakly on spectral type (cf. discussion in Hillenbrand et al. 1998). Accretion has the effect of 
filling in the line and the EW is then related to the mass accretion rate (Muzerolle et al. 1998). 
The observational values we will use are those obtained through 5 — 8A resolution spectroscopy 
by Hillenbrand (1997) for a large fraction of our optical sample (see Table 5 of Paper I). Figure 
12 (upper panel) shows the distribution of the Ca II EW both for our optical sample and for the 
subset of HRC detected stars. The detection fraction as a function of EW is shown in the lower 
panel. 

In order to make the distinction between accreting and non accreting stars more evident and to 
minimize the effects of observational uncertainties 7 , we divided our stars in two groups: stars whose 



7 Hillenbrand et al. (1998) quote EW errors of ±0.5. 
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Fig. 11. — Log(Lx/Lt,oi) vs. Log(M/M & ) in four age ranges, with symbols and lines same as 
Figure 10. Lower portion of plot (right of legend) indicates masses at which the structure of stars 
of Log(Age) are expected to be fully convective (leftward arrows) and fully radiative (rightward 
arrows), according to SDF evolution. 

Ca II line is seen in absorption, with EW(C& II) > 1, and which therefore are inferred to have low 
or no accretion, and stars that show the line in clear emission, with EW(C& II) < —1, and which are 
undergoing significant accretion. We exclude from our analysis stars with — 1 < EW (Ca II) < 1, 
about half the sample for which we have Ca II EW information (see Figure 12). As the majority 
of stars in the ONC appear to possess a circumstellar accretion disk, judging from both Ca II line 
and IR photometric indicators (Hillenbrand et al. 1998; Lada et al. 2000; Muench et al. 2001), we 
refer to stars in the two selected groups as having high and low circumstellar accretion and avoid 
using the terms "classical" or "weak-line" T Tauri stars. 

In Figure 13 we present X-ray luminosity functions for stars in the six lowest mass ranges 
defined in §2.2, computed separately for the two accretion-differentiated samples. Figure 14 shows 
median Log(Lx) vs. Log(M/M®) depicted in the same format as Figure 4 and separated into 
low- and high-accretion subsamples. For any mass bin with enough points, it is clear in either 
representation that stars with high accretion have lower median Lx than low-accretion stars. The 
difference of the medians is about an order of magnitude in the 0.25 — 2.0 M/M Q mass range. Two- 
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Fig. 12. — Upper panel: Distribution of Ca II (A = 8542^4) line equivalent width for the optical 
sample and the subset of stars detected in our HRC data. Lower panel: detection fraction. 



population statistical tests confirm the difference with high confidence (null-hypothesis probabilities 
less than 10~ 4 , for M = 0.25 — 0.50 M/M Q ). Similar plots for Lxj^boi confirm the same result, 
also with high confidence: low-accretion objects are more X-ray active that similar high-accretion 
objects. In order to check that the dependence of activity on accretion is not due to different 
mean ages/bolometric luminosities in the two subsets, we repeated this analysis in several mass 
and age slices, similar to those of Figure 10. When we did this, we invariably confirmed that the 
high- accretion stars are less X-ray active than the low-accretion ones. For any mass and age bin 
where the two subsamples are well represented, we obtain highly significant results (e.g., in the bin 
with M/Mq = 0.25 - 0.50 and Log(Age) = 6.0 - 6.5, null-hypothesis tests give better than 99.96% 
confidence that the two accretion categories differ). 

Differences in X-ray activity level aside, trends of mean activity with mass seem to hold 
preferentially for low accretion stars. The scatter in X-ray activity levels in any given mass range 
also appears lower for low-accretion stars than for the whole sample (cf. Figures 4 and 14). A 
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Fig. 13. — Optical sample XLFs for stars with high- and low-accretion (solid and dashed lines, 
respectively). Panels give different mass ranges, as indicated by legends. Legends also give numbers 
of detected (d) and undetected stars (u) used for XLFs of high- and low-accretion subsamples. 

relationship of activity with mass is not apparent for high-accretion stars, which also evidence 
larger spreads in their activity levels. 

3. The nature of unidentified sources 

Forty of our HRC sources remain unidentified with any optical/IR object. In this section we 
will explore the nature of these sources. 

As noted in Paper I (cf. §2.2), our X-ray detection method will yield about 10 spurious 
detections, and virtually all of these are expected to be unidentified (Paper I, §3.2) and have near- 
threshold SNRs. We therefore have ~ 30 (= 40 — 10) sources with no optical/IR identification 
("non-IDs") that are expected to be real and require some explanation 8 . Particularly puzzling is 



8 Of the 40 "non-IDs" 14 were previously detected with Chandra by Garmire et al. (2000) and/or Schulz et al. 
(2001); we are therefore confident that these 14 X-ray detections are not spurious. 
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Fig. 14. — Log(Lx) vs. Log(M / M Q ) for high-accretion (black) and low-accretion (gray) stars. 
Symbols have same meaning as in Figure 4, but with individual data points omitted for clarity. 
Because only three detections are found in the lowest mass bin of the high-accretion sample, we 
show an upper limit to the mean Lx in place of the median. 

the fact that most of them appear in the central 5' x 5' of the FOV (see Figure 15) where the 
-fT-band coverage is quite deep (limiting K ~ 17.5). A first clue to their nature comes from their 
spatial distribution and from the correlation with the Goldsmith, Bergin, & Lis (1997) 13 CO data 
on the molecular cloud total extinction, as shown in Figure 15. If we exclude four, likely-spurious 9 
sources with SNR < 4.8, the median V-band total extinction at the positions of the remaining 
sources is ~63 mag, much larger than that for all detected sources (~24 mag) in this central region. 
This suggests that the "non-IDs" as a group are high- intrinsic X-ray luminosity objects, which are 
either highly embedded in the cloud or which lie in the background. 

This conclusion is supported by our analysis of the hydrogen column density, Njj, measured 
directly from X-ray data (cf. Appendix A of Paper I). Since the "non-IDs" are largely low- 
SNR detections 9 , however, the hardness ratios derived from Advanced CCD Imaging Spectrometer 



9 We note that the 26 (=40 - 14, see footnote 8) unidentified detections that have not been confirmed have extremely 
low SNRs (median SNR = 4.94, vs. 12.45 for those previously confirmed by Chandra). 
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Fig. 15. — Spatial distribution of "non-IDs" (pluses) in the inner ~ 5' x 5' of the HRC field; squares 
indicate those previously detected by Chandra (Garmire et al. 2000; Schulz et al. 2001). Shaded 
contours trace molecular cloud total optical extinction (from Goldsmith, Bergin, & Lis 1997), with 
gray scale running from Ay < 10 (white) to Ay > 70 (dark grey) in steps of 10 magnitudes. 



(ACIS, Townsley et al. 2000) spectra are subject to large uncertainties. 10 Nevertheless, the fact 
that the low-energy portion of their spectra are extremely underpopulated with respect to higher 
energies is consistent with the hypothesis of high extinction: on average, "non-ID" hardness ratios 



The hardness-ratio approach sutlers not only from low source statistics but also from the fact that the low-energy 
hardness ratio saturates at high values of extinction (see Paper I, Appendix A) . 
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suggest optical extinction > 30 magnitudes. 

Having established that a majority of our unidentified sources are associated with highly ab- 
sorbed objects, we next consider whether or not the hypothesis that they are high-intrinsic lumi- 
nosity, embedded cluster members is compatible with the observational constraint that they are 
not detected in the deep K-band survey. We can derive an empirical relationship between stellar 
mass and intrinsic, non absorbed, i^T-band magnitude 11 by considering the two quantities for ONC 
members with low measured optical extinction (Ay < 1.0). We then derive, as a function of mass, 
the minimum iT-band extinction needed to make these sources fall below the JT-band sensitivity 
threshold of K = 17.5. Converting this minimum Ak to a minimum Njj = 9.3Ak ■ 2 • 10 21 cm~ 2 
and assuming a detection threshold of 10 photons, we then find the minimum Lx that stars emit- 
ting at a given temperature and similar to those in our low-extinction sample should have in order 
to be both detected in our HRC data and too absorbed to appear in the deep fT-band survey. 

Figure 16 shows the result of this exercise for three different assumed source temperatures 
within the range determined in Paper I. Any source that we detected without a i^-band counterpart 
requires an X-ray luminosity in excess of 10 30 3 erg-s -1 . The large scatter in inferred Lx values (for 
a given temperature) is due to the spread in the mass - K-band magnitude relation. We also show 
on the same plot the relation obtained in §2.2 between stellar mass and Lx for the optical sample. 
We conclude that "non-IDs" may well be embedded ONC members but must have masses larger 
than ~ 1 Mq, i.e. lower mass ONC members of the same activity level as our optical sample, if 
detected in X rays, must appear in the K-band survey and hence cannot be "non-IDs." 

Thus we have established that the "non-IDs" are indeed affected by high extinction and that 
their HRC count rates and ET-band magnitude lower limits are compatible with their being highly 
absorbed "high" mass PMS stars similar to those of our optical sample. In order to confirm their 
ONC membership, we must also exclude the possibility of their being background extragalactic 
objects. 

We can estimate the number of expected extragalactic X-ray sources by combining information 
on the expected sky density of such objects, as a function of unabsorbed flux (Giacconi et al. 
2001; Tozzi et al. 2001), with our knowledge of the total extinction due to the Orion Molecular 
Cloud (Goldsmith, Bergin, & Lis 1997). For this exercise we make the following assumptions: 1) 
Extragalactic sources have power law spectra of photon index T between 1.3 and 1.5, compatible 
with values derived by Giacconi et al. (2001); 2) The sky density p[F] of extragalactic objects as 
a function of unabsorbed flux follows the descriptions of Giacconi et al. (2001) and Tozzi et al. 
(2001) 12 ; 3) Observed photon flux is determined by intrinsic flux and an absorption proportional to 



11 There is significant scatter in K for this relation, probably due both to variable IR excesses induced by the 
presence of circumstellar disks and to a spread in stellar ages. 

12 They give two relations, p s {F s ) and ph(Fh), for the sky density of objects with unabsorbed flux larger then F, 
measured in soft (0.5-2.0 keV) and hard (2.0-10.0 keV) spectral bands, respectively. 
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Fig. 16. — Log-log plot vs. mass of inferred minimum X-ray luminosity for objects too faint to 
appear in the K-band survey of Hillenbrand & Carpenter (2000), yet X-ray luminous enough for 
Chandra detection. Triangles give hypothetical lower limits for stars in our optical sample with 
Ay < 1 (see text), computed for different values of X-ray temperature as indicated by legend. For 
reference, we also show the characterization of the XLF as presented in Figure 4. 



the molecular cloud total optical absorption (Njj = 2 • lO 21 ^?'), as estimated from radio data by 
Goldsmith, Bergin, & Lis (1997); we use CF(Nh) to denote the unabsorbed flux to HRC count rate 
conversion factor computed with PIMMS, which depends on Nh, the assumed source spectrum and 
the spectral band under consideration; and 4) We can detect sources above count rate thresholds, 
CRth, computed as the photon threshold (viz. 6, 8 or 10) divided by exposure time. With these 
assumptions, the surface density of extragalactic sources seen through a given hydrogen column 
density becomes p[CR t h • CF(Nh)] and we can compute the expected number of extragalactic 
sources detected in an area S as: 



/ p[CR th ■ CF{N H )\ ■ da 
Js 



(1) 
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where da is the infinitesimal area element and Nh is a function of position on the sky. Computing 
the integral (1) for the given values of T, CRth and energy band (and therefore /9[-F]), we find 
expected numbers of detectable extragalactic objects in the inner 5' x 5' of our FOV ranging from 
0.6 to 1.2. An extragalactic nature for most of our unidentified sources is therefore excluded, and 
most "non-IDs" can indeed be highly embedded ONC members never observed with optical or 
near-IR instruments. 

We can get a rough estimate for the size of this hidden population if we assume that coun- 
terparts of the ~20 unidentified field-center sources (expected to be real) have masses larger than 
~ 1 M Q (cf. Figure 16) and that the ratio of the number of stars with M > 1.0 M to stars with 
M < 1.0 Mq is the same as for our optical sample, namely, ~25% for the field center and ~16% 
for the entire FOV. We obtain the interesting result that the undiscovered population numbers be- 
tween 80 and 125 stars, i.e., about half of the optically- well characterized members 13 in the same 
region and about a quarter of the optical objects studied by Hillenbrand (1997). For comparison, 
there are ~ 600 stars seen in the infrared with K < 14.0, implying masses that roughly correspond 
to the mass range spanned by our optical sample. We conclude, therefore, that the central region 
of the ONC contains ~100 additional stars that remain hidden from the deepest -ftT-band surveys 
conducted to date. 



4. Discussion and Summary 

Using Chandra X-ray data and recent optical/IR data from the literature, both presented in 
Paper I, we have studied relationships among X-ray activity indices and various stellar parameters 
for an optically-selected sample of ONC members. Our ultimate goal is to explain on physical 
grounds the mechanisms responsible for X-ray emission and magnetic activity in PMS stars. 

We do not find any correlation of either Lx or Lx/L^i with stellar rotational period. This 
result has sometimes been regarded as evidence against the a - u> dynamo explanation of activity 
for PMS stars. We note however that: 1) Although the sample of stars for which we have a 
measure of rotational period is larger than ever considered in similar past studies, it is still biased 
toward more active stars, as judged from distributions of X-ray activity indicators; 2) Many of the 
sources with known rotational properties may thus have saturated activity levels and the intensity 
of their X-ray emission may therefore be insensitive to the rotational period; 3) 90% of these 
rotational periods are shorter than 10 days; 4) As discussed by Flaccomio (2002), periods up to 
~ 16 days might correspond to saturated activity for ~1 Myr-old stars, according to the consolidated 
activity - Rossby number relation. We therefore do not regard the present lack of evidence of a 
rotation-activity relationship as conclusive for the ONC and look forward the realization of more 
complete rotational period databases. We also note that even if the measured rotational periods 



About 200 stars are placed in the HR diagram, and ~140 are in our optical sample. 
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were representative of the whole ONC population (as recently indicated by Rhode et al. 2001), a 
lack of correlation does not necessarily require a new explanation of activity 

As previously noted for other star forming regions (e.g., Feigelson et al. 1993), we find a 
significant correlation of Lx with mass for M < 3.0 Mq. The logarithm of the median Lx rises 
quite linearly on a logarithmic mass scale from Log{M / Mq) ~ —0.9 to ~ 0.4, but for higher masses, 
drops sharply before rising again above Log(M/M & ) > 1.0. The drop in activity at M ~ 3.0 Mq is 
even more dramatic in Lx / L^ i (more than two orders of magnitude) and may be further enhanced 
if proper account for multiplicity were taken. Supported by recent PMS evolutionary models (SDF) 
we interpret this drop as the effect of the disappearance of a convective envelope for stars more 
massive than ~ 3 Mq. In general we note a convincing correlation between the presence of a 
convective layer and activity for M < 10 Mq. 

For 0.5 < M/Mq < 3.0, the median Log{Lx / L^i) is observed to be quite stable at about -3.5, 
lower than (though close to) the saturation level. At lower masses however we observe a decrease 
of the median Lx/L^i with decreasing mass, as it falls by a factor of ~ 2.8 (0.45 dex) between 
M = 0.5 — 1.0 and M = 0.1 Mq. We are unsure of the origin of this mass dependence but note a 
rough coincidence between the turning point in the Lx/L^i - mass relation, at ~ 1 M & , and the 
mass at which 1 Myr stars cease to be fully convective (cf. Figure 11). We suggest this may be due 
to a transition from an a—co dynamo to a fully-convective dynamo. It may be even more relevant to 
note that the median value of the Ca II line equivalent width, which can be used as a proxy of disk 
accretion (cf. §2.5), is significantly lower for stars with masses below 0.5M Q (median EW=-0.15) 
than for stars with mass between 0.5 and 3.OM (median £W=0.94), possibly indicating a larger 
fraction of accreting stars at low masses 14 . We therefore propose (cf. Flaccomio 2002) that the 
increase of Lx/L^i at low masses might be the effect of a decrease, with increasing stellar mass, 
of the fraction of stars that are accreting and/or are surrounded by disks. We indeed found such 
stars to have significantly lower activity levels (§2.5). 

Although we do not detect any securely identified brown dwarfs, we are able to estimate the 
mean X-ray luminosity for the 9 spectroscopically confirmed BDs of the BD bin in Figure 7. The 
mean Lx of ~ 10 28 ' 5 erg-s -1 appears to be on the same Lx vs. mass relationship we observe for 
higher mass stars. This is not particularly surprising because at this evolutionary stage BDs do 
not fundamentally differ from low mass fully convective PMS stars. 

We have evidence of the presence of a population of very weak emitters at the low mass 
end. In the 0.25 — 0.50 Mq mass bin, for example, we observe that the mean X-ray luminosity 
of stars individually not detected in the X-ray data (Paper I) is > 1.5 orders of magnitude lower 
with respect to the mean of the whole sample, indicating either a very wide or a bimodal X-ray 
luminosity distribution. At least part of this dispersion can be explained by a dependence of 



There is indeed growing evidence that disk lifetime is longer for lower mass stars. Two examples of such a result, 
based on detection of disks in the IR, can be found in Hillenbrand et al. (1998) for the ONC and Haisch et al. (2001) 
for IC 348. 



-24- 



activity on circumstellar accretion rate. By discriminating stars with high and low circumstellar 
accretion on the basis of the Ca II (A = 8542A) equivalent width, we prove with high statistical 
significance that the X-ray activity of the two groups differs by as much as an order of magnitude. 
The relationships of median Lx and Lxl^boi with mass appear to hold exclusively for the low 
accretion sample. We also note that the scatter in activity at a given mass appears to be larger 
for the high accretion sample than for the low accretion sample. Although the origin of these 
differences is presently unknown, we put forward three possible scenarios: 

1) Accreting stars are rotating slowly due to disk breaking and their dynamo efficiency depends 
on rotation, while those accreting less have broken free of their disks to spin-up and thus saturate 
their activity. This hypothesis seems to be contradicted by two facts: low and high accretion stars 
have statistically indistinguishable rotational periods (cf. Figure 1) and, as reported in §2.1, we do 
not observe a correlation between activity and P ro t- However, the sample of stars for which we have 
rotational information is not complete and may be subject to biases (see §2.1), so that we cannot 
rule out this possibility until a more representative sample of rotational periods becomes available. 

2) Accretion and/or the presence of a disk, and/or outflows influences coronal geometry, for 
example, by decreasing the fraction of the stellar surface available for the closed magnetic structures 
from which X-ray emission is thought to emanate. Other scenarios of altered geometries have been 
proposed; e.g., Montmerle et al. (2000) argue that coronal structures might extend to the inner part 
of disk. Because of the inhomogeneous and time variable nature of accretion, variability studies 
and simultaneous X-ray/optical observations could help clarify this matter. 

3) Accreting stars have higher X-ray extinction than assumed, so that the difference in inferred 
Lx and Lx/L^i is only apparent. Although in Paper I our ACIS spectra analysis was not able 
to exclude this possibility with high confidence, due mainly to the low photon statistics of sources 
associated with high accretion stars, we also found no supporting evidence for this scenario. More- 
over, an empirical determination of the Nh vs. Ay relation for the p Ophiuchi population (cf. 
Imanishi et al. 2001) also does not support this hypothesis. The reason for the departure from the 
assumed Ay-Nn relation could be that the gas to dust ratio may differ for accreting stars with 
respect to the average interstellar value, resulting in an underestimation of the X-ray absorption. 
One could imagine a scenario in which accretion gas columns that cross the line of sight would 
obscure X-rays and let optical/IR radiation through. The most obvious test for this hypothesis 
is a better determination of hydrogen absorbing columns through deeper (i.e., longer exposures) 
medium spectral resolution X-ray observations. As for the previous hypothesis, variability studies 
might also give useful clues. 

Finally we have investigated the nature of sources detected in our HRC X-ray data and not 
identified with any optical/IR object. On the basis of their X-ray spectra and location in the cloud 
we conclude that, as a class, they are subject to high extinction. We exclude an extragalactic nature 
for most of these objects by considering the sky density of such objects as a function of X-ray flux 
and of the molecular cloud total extinction and conclude that most unidentified X-ray sources are 
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likely associated with deeply embedded, high intrinsic X-ray luminosity, ONC members. Given the 
relationship between Lx and stellar mass we infer that these stars are most likely fairly massive 
(M > 1.0 Mq), thus indicating the likely presence of a significant number of lower mass, as yet 
undetected ONC members. 

In light of these results, we propose this scenario: 

1) Activity in ~ 1 Myr old low mass PMS stars is ultimately due to a dynamo mechanism 
that requires a convective layer in order to function. Rotation may or may not be a fundamental 
ingredient of this dynamo. 

2) The same mechanism responsible for X-ray emission of low mass stars is probably also at 
work in brown dwarfs. 

3) Many of our stars appear to be saturated. The fraction of saturated stars, or alternatively the 
saturation level, depends on stellar mass. At low masses we might be seeing the effect of decreasing 
accretion/disk fraction with increasing mass, or the transition between two kinds of dynamos: one 
that functions on fully convective stars and brown dwarfs and in which the saturation level depends 
on mass, and another at work on partially convective - partially radiative stars, characterized by a 
constant saturation level. 

4) Highly accreting stars are seen to be less active with respect to low accretion counterparts, 
possibly because of (a) different rotational properties, (b) disk/accretion induced modifications of 
the coronal magnetic field geometry, or (c) anomalous interstellar X-ray absorption due to the 
presence of circumstellar material. 

Additional deep X-ray and IR observations are needed to advance our understanding and 
enable us to identify which scenario(s) Nature prefers. 
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A. "Composite" source analysis 

Here we describe a "composite" source method we employed to study the mean X-ray properties 
of groups of ONC members, some or all of which were individually undetected in our HRC data. 
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The method was applied in §2.2 to undetected low mass stars, and in §2.3, to objects of substellar 
mass. 

The basic steps of our method are to: 1) define a suitable sample of stellar or substellar 
objects; 2) determine the number of photons detected in a circular region at the X-ray position 
corresponding to each optical object; 3) correct the photon numbers for background contributions 
and the fact that a (predictable) fraction of the photons from a point-like source will fall outside 
the extraction circle; 4) sum the contributions of single objects and derive a signal-to-noise ratio 
(SNR) for the "composite" source; 5) divide the "composite" source net counts (or upper limit, 
for SNRs below a set threshold) by the cumulative exposure time to obtain a mean count rate for 
objects in the sample; and finally, 6) convert the mean count rate (or upper limit) to a "composite" 
Lx value. 

The critical points in these six steps for estimating Lx values for the composite sources are as 
follows: 

1) The selection of composite source samples was based on stellar (or substellar) mass as given 
by our optical sample. We also added three more requirements: a) that the objects lay in the inner 
5' of the HRC field of view, b) that the closest detected X-ray source (other than the object in 
question, if detected) was farther away than 5" and c) that the measured optical extinction Ay was 
less than 2.0. Constraint (a) ensures that the PSF at the position of all the objects in the sample 
contains more than ~ 60% of the photons falling within a circle of 2" radius, while (b) isolates the 
source to avoid contamination with photons from neighboring sources. Constraint (c) reduces the 
uncertainty in the conversion between composite source count rates and X-ray luminosities. 

2) Individual source photons were extracted from a 2" circle centered on the optical star po- 
sition. Figure 17 shows an example of the extraction regions for undetected members in the 
0.16 — 0.25 M/Mq mass range as well as a graphical representation of the composite source resulting 
from the summation process. 

3) The background contribution in each of the extraction circles was estimated from the values 
of the source- free background map produced by pwdetect (cf. Paper I). The fraction of counts 
from the object falling within this extraction circle ranges from ~ 60% to 95% and was estimated 
from the images of calibrated PSFs obtained for the position of each source using the CIAO task 
mkpsf and the Fl standard PSF library file (Karovska et al. 2001). Sources were assumed for this 
purpose to be monochromatic with an energy of 2.0 keV. A visual comparison of the expected and 
observed shape of the PSF for bright sources in our present sample indicates a good qualitative 
match. 
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4) From counting statistics we estimated the probability that the measured total number of 
extracted counts is the result of random fluctuation in the background level. In those cases in 
which such probability was larger than 0.135% (3<r), we computed upper limits to the composite 
source counts using the same confidence threshold. 

5) The cumulative exposure time used to compute mean count rates for the objects in each mass 
segregated sample was obtained by summing the values of the exposure map (cf. Paper I) at the 
position of the objects. 

6) We converted mean count rates to X-ray luminosities in a manner consistent with Paper I, but 
such conversion depends critically on extinction. While we can derive conversion factors for each of 
the component objects, a global conversion factor for a composite source cannot be readily defined. 
If all the objects in the sample had the same intrinsic X-ray luminosity the conversion factor for 
the average source would be the average of the single conversion factors. In order to reduce the 
range of individual conversion factors we have limited our object sample to low optical extinction 
{Ay < 2.0). However some uncertainty remains: as shown in Figure 8 of Paper I, this range of 
Ay corresponds to a factor of about 3 (~ 0.5 dex) variation in the conversion factor at our chosen 
source temperature (2.16 keV). We therefore report two values of luminosity for each of our average 
sources: one computed with the mean of the conversion factors, the other with the median. 

The results of applying this technique are discussed in §2.2 and §2.3 and are shown in Figure 

7. 
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Table 1. Detections and Upper Limits for mass/age ranges 



Log(Ag 


e ) [y r -] 


0.1-0.16 


0.16-0.25 


Mass [Mq] 
0.25-0.50 


0.5-1.0 


1.0-2.0 


2.0-3.0 


3.0-10.0 


4.5-5.5 


Det. 


1 


14 


13 


3 





1 





11 


Upp. 


6 


12 


6 














5.5-6.0 


Det. 


3 


7 


27 


19 


1 


2 


3 


11 


Upp. 


13 


16 


18 


6 








1 


6.0-6.5 


Det. 


3 


21 


95 


31 


12 


18 


2 


11 


Upp. 


8 


43 


80 


10 


3 


1 


3 


6.5-7.5 


Det. 


3 


7 


7 


10 


20 


7 





11 


Upp. 


24 


39 


27 


5 


4 


7 






